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Abstract 

We prove an inequality for the spectral radius of products of non-negative matrices 
conjectured by X. Zhan. We show that for all n x n non-negative matrices A and 
B, p(A o B) < p((A o A)(B o B)) 1 / 2 < p(AB), where o represents the Hadamard 
product. 
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We denote by Aij the entry of a matrix A in position and by o the 

Hadamard product, that is, the entrywise product (AoB)^ = AijBij. We shall 
be considered with entrywise positive and entrywise non-negative matrices, 
which we'll simply call positive and non-negative matrices, and denote by 
A > and A > 0, respectively. 

The spectral radius p(A) of a matrix A is the largest modulus of the eigenvalues 
of A. For non-normal matrices, the spectrum does not build norms, for the 
simple reason that a matrix can be unbounded even when all its eigenvalues are 

(oi\ 

zero. The Jordan block serves as the prime example here. In particular, 

the spectral radius is not a matrix norm. This shows up for example in the 
fact that the spectral radius is not submultiplicative: p(AB) < p(A)p(B) does 
not hold in general, not even when restricting to non-negative matrices. A 
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counterexample is given by the pair [3] 





A 



On the other hand, for non-negative A and B, the spectral radius is submulti- 
plicative w.r.t. the Hadamard product: p(AoB) < p(A)p(B) ([3], observation 
5.7.4). 

In the light of the few relevant statements that can be made about the spectral 
radius, X. Zhan made the remarkable discovery that for pairs of square non- 
negative matrices, the spectral radius of the Hadamard product is always 
bounded above by the spectral radius of the (conventional) matrix product. 
This was presented in the talk |lj and posed as a conjecture. 

The purpose of this note is to prove this conjecture. In fact, we will prove a 
little more: 

Theorem 1 For n x n non-negative matrices A and B, 

p{A oB)< p{{A oA){Bo B)) 1/2 < p{AB). 



Before going into the proof, we note that there is no reasonable lower bound 
on p(AoB) in terms of p(AB). The pair A and B given above also shows that 
p(A o B) can be zero while p{AB) is not. In fact, this can happen even when 
A and B both have non-zero spectral radius. Take 



A 



1 + xy x 

Then A, B and AoB all have spectral radius 1, while AB = \ ] . Since 

y i 

det(AB) = 1 and Tr(AB) = 2 + xy, the spectral radius of AB is unbounded. 




Let us now present the promised proof of Theorem 1. The proof is based on 
the following representation of the spectral radius of positive matrices: 

Lemma 1 For A e M n such that A > 0, p(A) = lim m ^ 00 (Tr A m )^ m . 

Proof. By a basic theorem in Perron-Frobenius theory (see e.g. Theorem 8.2.8 
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in [2]), for any A > we have 



( 



.4 



) 



m 



lim 



T 

= xy , 
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where x and y are certain positive vectors for which x T y = 1. Taking the trace 
yields 



and the lemma follows immediately □ 

The lemma is not generally true for non-negative matrices, because the limit 
need not even exist. The simplest example is the matrix 



whose spectral radius is 1. While its even powers are equal to the identity 
matrix and have trace equal to 2, its odd powers are equal to A itself and 
have trace equal to 0. More generally, every permutation matrix shows this 
kind of behaviour, with a periodicity given by the least common multiple of all 
the cycle lengths of the corresponding permutation. In these cases, replacing 
the limit by the limit superior would salvage the lemma, but we do not know 
whether this will be true in general. On the other hand, the statement of the 
lemma does still hold unmodified for primitive non-negative matrices. 

However, for our purposes we do not need all these generalisations. In fact, 
we will only explicitly prove p(A o B) < p(AB) for positive A and B. Because 
of the continuity of the spectral radius, the Theorem then follows for non- 
negative A and B as well. 

An obvious consequence of the lemma is that for A > 0, and any positive 
integer k, p(A k ) = p(A) k . By continuity of the spectral radius, this is also true 
for A>0. 

Proof of Theorem 1 for positive A and B. We will prove that the following 
inequalities hold for any positive integer k: 





Tr((A o B) 2k ) < Tr((y4 o A){B o B)f < Tr((AB) 2k ). 



By taking the (2/c)th root and taking the limit k 
statement of the theorem, by the lemma. 



oo, this implies the 



The left-hand side Tr((A o B) 2k ) can be written as a 2A;-fold sum: 
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Tr((A o B) 2k ) — ^2 (A ili2 Bi li2 ) (A i2is B i2i3 ) . . . (A i2kil B i2kil ) 

— ( y 4jli2-^«2«3 • • • ^■«2fe-li2fc-^«2fe«l) 

H,«2vi*2k 

X (-Sil«2^-«2«3 • • • -^«2fc-l*2fc^2fcn) • 

Note that the alternation of the A and B factors in the last line is intentional. 



The last expression can be seen as an inner product between two vectors 
in R™ , one with entries A ili2 Bi 2 i 3 . . . Ai 2k _ li2k Bi 2kil and the other with en- 
tries B ili2 Ai 2 i 3 . . . Bi 2k _ x i 2k Ai 2k i x . One sees that these two vectors are the same 
up to a permutation of the entries (as can be seen by performing a cyclic 
permutation on the indices ix, ?2, • • • , izk)- Thus, in particular, both vectors 
have the same Euclidean norm. Applying the Cauchy-Schwarz inequality, 



(x,y) < J{x,x) (y,y), then gives 



Tr((A o B) 2k ) < £ (Ai li2 B i2i3 . . . A i2k _ lt2k B t2kll 

(^Ai 1 i 2 Bi 2 i. i . . . Ai 2k _ 1 i 2k Bi 2k i 1 



ll,l2i".,»2fc 
X 



Noting that the last line can be succinctly written as Tr((A o A)(B o B)) k 
proves the first inequality of the theorem. 

Now consider Tr((AB) 2k ). This can be written 4fc-fold summation: 



Tr((AB) 2fe ) — £ (Aii 2 -^i2«3 • • • Ai2k-ii2k-Bi 2k ji) 

n.J2>--vJ2fc 

x ^hh^hh ■ ■ ■ ^hh-ihh^hhii) ■ 

The crucial observation is that if we take all terms of this summation for which 
h = hi ^2 = Hi - ■ -iik = 3iki then we obtain the right-hand side of p]). The 
terms we left out are of course all positive, because all matrix elements of A 
and B are positive. Therefore, we find that Tr((AB) 2k ) is an upper bound on 
(Tjrj, and thus on Tr((v4 o B) 2k ). This proves the second inequality. □ 
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